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We analyze the random Euclidean bipartite matching problem on the hypertorus in d dimensions 
with quadratic cost and we derive the two-point correlation function for the optimal matching, using 
a proper ansatz introduced by Caracciolo et al. [I] to evaluate the average optimal matching cost. 
We consider both the grid-Poisson (gP) matching problem and the Poisson-Poisson (PP) matching 
problem. We also show that the correlation function is strictly related to the Green’s function of 
the Laplace operator on the analyzed manifold. 


I. THE MATCHING PROBLEM: AN 
INTRODUCTION 

The assignment problem, or bipartite matching prob¬ 
lem, is a classical combinatorial optimization problem in 
which two sets of N elements, B = {bi}i=i,...,jv and 
TZ := {ri}i=i,...,Ar are considered. An assignment is an 
element ir G Sn oi the set iSjv of permutations of N el¬ 
ements, such that Vi Moreover, a cost function 

w: TZx B ^ R+, (r^, bj) >->• Wij, is given, from which we 
can define the total cost of a certain assignment tt as 

1 ^ 

En[tt;w]:=—'^W i„(^i). (1) 

i=l 

In the bipartite matching problem we want to find the 
permutation tt* that minimizes the previous quantity for 
a given function w. From a computational point of view, 
the problem belongs to the P computational complex¬ 
ity class and it can be efficiently solved using fast algo¬ 
rithms [m]. If random instances are considered, i.e. Wij 
are random quantities, we are usually interested on the 
average optimal cost, 

En ■= EN['7r*-,w] ( 2 ) 

where we denoted by • the expectation over all the possi¬ 
ble instances w. If the values {wij} are independent and 
identically distributed random variables, the problem is 
usually called random assignment problem: in this case, 
the average optimal cost and its properties in the large 
N limit were investigated both with statistical physics 
techniques [5] and probability arguments [B]. 

In a more complicated variation of the random assign¬ 
ment problem, the so called Euclidean bipartite match¬ 
ing problem (Ebmp), the two sets TZ and B are in one- 
to-one correspondence with uniformly generated random 


points on the unit hypercube Fla •= [0; C IR'^, whilst 
the weight Wij is a function of the Euclidean distance 
||ri — bjll (for simplicity, we identify the elements b^, Vj 
with the corresponding geometric points in Fid): in this 
case correlations between different values {wij} appear 
and a proper mathematical treatment is more compli¬ 
cated. In the following we will consider weight functions 
in the form 


Wij := \\r, - hj\f, pGR+. (3) 

Mezard and Parisi [7] considered the previous problem 
on the hypercube Fid for any value of d in the large N 
limit, assuming that correlations can be treated as per¬ 
turbations to the purely random case and evaluating ap¬ 
proximately the average optimal cost through replica ar¬ 
guments. An exact solution to the problem for the d = 1 
and p > 1 case is provided in mis], where the average op¬ 
timal cost and correlation functions are computed; more¬ 
over, the correspondence between the matching problem 
and the Brownian bridge process on the line and the cir¬ 
cle is proved. 

Denoting by tt* the optimal permutation for a given 
instance, we introduce the optimal matching ray: 

m(ri) := b^.(i) - r,, i = l,...N. (4) 

The optimal cost is 

<)[m;{7^,R}]:=iEtlllm(u)|^, (5) 

<)(d) :=<^[m;{7^,R}] (6) 

where d is the dimensionality of the Euclidean space 
and the average • is performed over the positions of the 
points. 

The scaling properties of the optimal matching ray, and 
therefore of the optimal cost, are known to the literature 
for p > 1 HHIl], being 


[A for d = 1, 

ford =2, 
\lkf for d > 3. 
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||m(x)|| 


(7) 
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Figure 1. gP Euclidean bipartite matching with N = 225 and 
p = 2 on the torus. 


In the present paper we are interested in the correlation 
function of the optimal matching ray m in the large N 
limit; we will assume periodic boundary conditions on the 
hypercube fid, i.e., we will consider the problem on the 
flat hypertorus in d dimensions. We will analyze both 
the case in which two sets of random points are consid¬ 
ered and the case in which one set of points is supposed 
fixed on a regular hypercubic lattice, whilst the second 
set is obtained from a Poisson process (see, e.g., fig. |l] 
for a pictorial representation of a realization of the two 
dimensional problem). To introduce our results, in Sec¬ 
tion [IT] we will review the Monge-Kantorovic formulation 
of the optimal transport problem, from which a suitable 
general ansatz, already used by Caracciolo et al. m , is 
derived for the expression of the optimal matching ray in 
the continuum limit for the p = 2 case. Using this work¬ 
ing ansatz, in Section |II B[ we will consider the Ebmp 
with quadratic cost oiiT“, and we will give evidences 
that, in the large N limit, the correlation function 


C'd(x) 


m(ri) • m(rj) 


-rj=x 


( 8 ) 


is related to the Green’s function of the Laplacian opera¬ 
tor on T'’*. We will consider also the correlation function 
for the normalized optimal matching ray 


cr(r,) 


m(ri) 

llm(ri)|| ■ 


(9) 


Finally, we will give numerical evidences that the func¬ 
tional forms of the correlation functions obtained for 
p = 2 in the two dimensional case are in good agree¬ 
ment with the numerical results for p = I and p = 3 in 
the same dimension. 

To our knowledge, these results are new to the liter¬ 
ature, where only the d = 1 case is evaluated explicitly 
BM- The present work can be seen as a natural ex¬ 
pansion and completion of a previous work of Caracciolo 
et al. [T]. 


II. THE MONGE-KANTOROVIC MASS 
TRANSFER PROBLEM AND CORRELATION 
FUNCTIONS 


A. The Monge-Kantorovic problem 

The Monge-Kantorovic transfer problem is a well stud¬ 
ied problem in Measure Theory [12], in the context of 
Transportation theory. Let us suppose that we are given 
two nonnegative smooth density measures pi and p 2 on 
a compact subset O C E,'’*, Pi(x) d'^ x = P 2 (x) x. 

Suppose also that a work function w: fl x fl —i' E+ is 
given. We want to find a one-to-one mapping M: fl —>■ fl, 
called optimal transport map, such that the following con¬ 
ditions are satisfied: 

1. M S T, T set of suitable transport maps, defined 
as follows 

r := |M: fl U|/^pi(x) d'^a: = P 2 (x) d'^ a: 

VA C fl measurable} . (10) 


2. M minimizes a certain cost functional 

E[M.-,w\:= w(x, A4(x))pi(x)d‘^a:, (11) 

where w: fl x fl R+ is a transportation cost 
function, i.e., if[M; w] = minj^^q- E[M.-^w\. 

Interestingly, it can be proven that, if w(x, y) = ||x — y||^ 
with p e (l,-|-oo) the condition A4 e T can be re¬ 
expressed as a change-of-variable formula [T3| 

Pi(x) = p 2 (M(x)) det Jav4(x), Vx e U, (12) 

where Jaw(x) is the Jacobian matrix of Ad., (Jaw(x))y = 
^l^(x). Moreover, if p = 2, it can be proved [T3| that 
the optimal transport map can be expressed as a gradient 
of a scalar potential p, i.e., 

M(x) = Vp(x). (13) 


It follows that the solution M of the transport problem 
has to be identified among the vector fields in the form 
(13), with if satisfying the following relation: 


pi(x) = p 2 (Vp(x))detHessp(x), Vx G fl. (14) 


In the previous expression Hessp(x) is the Hessian ma- 
trix, (Hessp(x))„. = gw^(x). The nonlinear equation 


(14) is known to the literature as Monge-Ampere egua- 


tion. 


Consider now the Monge-Kantorovic problem in fl C 
X = Jq P 2 {x) d"^ a; = 1, with work function 
w{x,y) = |jx — y||^; suppose moreover that 


pi(x) = 1-1- Spi{x) and p 2 {x) = 1 -|- Sp 2 {x), (15) 
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where 


1(5/91 (x)I < 1 and |5p2(x)| <1 Vx e (16) 

We expect that M(x) = x4-m(x), ||m(x)|| 1 Vx G f2: 

in the first order approximation, det Jm(x) « 1 + V • 
m(x), so we have that 

V • m(x) « pi(x) - p 2 (x) =: p(x); (17) 

in particular, using the fact that m = W(p, we have that 
in the limit of our approximation the Poisson equation 
holds: 


A(/)(x) =p(x). (18) 


Note that in this case the total cost of the transport is 
given by 

if[M,||.f] = ^||V(/)(x)f d'^cr. (19) 


In the following we will consider the unit hypercube, 
n = with periodic boundary conditions, i.e., we 
will work on the flat hypertorus := Being 

p(x)d^x = 0, Eq. (18) has a unique solution on the 
compact manifold given by 


</> W = /jd p{y)Gd{y, x) d”^ y, 


( 20 ) 


where Gd is the Green’s function for the Laplace operator 
A on T'’* defined by the relation 


limit, the Ebmp between two sets of points on the hy¬ 
pertorus appears as a transport problem between two 
atomic measures that can be assumed as almost uniform 
measures on the domain of interest. This naive approach 
is justified a posteriori by the excellent agreement be¬ 
tween theoretical predictions and numerical results. In 
particular, let us denote by B := C Bid and 

TZ := C fid two sets of points in fid, each 

set of cardinality N. The optimal cost of the matching 
on the flat hypertorus T^^ with quadratic cost is therefore 
given, in the notation above, by 

1 ^ 

= — ^ ||m(r,)f. (26) 

' i=l 

In the previous formula, m is the (geodesic) optimal 
matching ray on T”^. Let us now introduce two atomic 
measure densities: 

Pn{y^) = j^Y.^=i5G){x-Yi), (27a) 

Pb(x) = (27b) 

In the following, we will assume, as a working ansatz, 
that, in the large N limit, m —>■ m(x) = V(()(x) in such a 
way that M(x) = x-|-m(x) is an optimal transport map 
between the two almost uniform measures above. Under 
the hypothesis that at least one set of points is randomly 
generated, we introduce the following correlation func¬ 
tion 


AyGd(x,y) =5(‘^)(x-y)-l, (21) 


Gd(x,y) := V(()(x) • V0(y) (28) 


the solution of which can be written as 


Gd(x,y) = Gd(x-y) 


E 

nGZ‘^\{0} 


g27rm.(x-y) 

47r2||nf 


( 22 ) 


In this linear approximation, the transport cost is given 

by 


E[M, ll-f] = - p{x)Gd{x,y)p{y) d‘^ y X 
|p(n)P 


(23) 

(24) 


where 


and, by using Eq. (20), we obtain 


Gd(x,y) = Gd(x-y) = 

= Jf ^ 2 :Gd(z - x) • VwGd(’w - y)p(z)p(w) d'^ z d'’* w, 

(29) 


where we denoted by 

p(x) := pTi{x) - pb(x) 

= (x-r,)-5('') (x-b,) 

2=1 


(30) 


/o(n) := fj^p(x)e ^^-^d'^a;. (25) 

B. Correlation functions for the EBMP on the 
hypertorus 

In the previous section we introduced the Monge- 
Kantorovic transport problem and we obtained also a set 
of simple results for the transport problem with quadratic 
cost between two almost uniform measures, through a 
proper linearization. We want to extract useful informa¬ 
tion about the discrete combinatorial problem from the 
continuum problem using the fact that, in the large IV 


and the average • is intended over all possible instances. 
Observe now that Gd(0) is, in the large N limit, the av¬ 
erage optimal cost for the Euclidean bipartite matching 
problem; using this simple correspondence, Caracciolo 
et al. [1] derived the correct scaling of the optimal cost 
and, through a proper regularization procedure, the fi¬ 
nite size corrections to the average optimal cost for any 
dimension. In the following we will consider the complete 
correlation function Gd(x) in any dimension and we will 
derive it using the same ansatz successfully adopted by 
Caracciolo et al. [T] to obtain the scaling of the average 
optimal cost. 

We will distinguish two different cases. 
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Poisson Poisson Euclidean matching problem: 

In the Poisson-Poisson (PP) Euclidean matching 
problem both the points of TZ and the points of 
B are random points uniformly distributed within 
rid- In this case we obtain 


PWp(y) = ^ 5(^)(x-y)-I 


(31) 


and therefore the correlation function is 


C'd(x-y) = - —Gd(x-y). 


(32) 


As anticipated, the average optimal cost is given by 


in which the correlation between the values normalized 
transport field 


cr(x) := 


m(x) 

|m(x)| 


Vx<?i'(x) 

|Vx0(x)|| 


(38) 


in different positions is evaluated. Note that cr lives on 
the d-dimensional unit sphere. To compute the correla¬ 
tion function (37) for the normalized field in the PP case, 


we assume a Gaussian behavior for the joint probability 
distribution of two values of the optimal transport field, 
and therefore we have 


:= {TZ,B}] = 0^(0). (33) 

Grid—Poisson Euclidean matching problem: In 

the grid~Poisson (gP) Euclidean matching prob¬ 
lem we suppose that N = L'^ for some natural 
number L S IN and that one set of points, 
e.g. the set TZ = {ri}j=i,...^Ar, is fixed on the 
vertices of an hypercubic lattice, in such a way 
that TZ = { ^|fc G (0, n IN'^}, whilst the set 
B = {bi}j^i_,.._jv C fid is obtained as before 
considering randomly generated points in Vld- We 
have 

p(x)p(y)=AdW(x-y) + ^ 

+ l^E.,<5(‘')(x-G)<5W(y-r,) (34) 

(x-r,) + dW (y-r,)]. 


cd(x - y) = 



mi ■ m2 

llmillllmall 


e 


(mi m 2 ).S-i(x,y)-(™J) 
^27r-\/det 

(39) 


where S(x, y) is the covariance matrix. 


f m(x) • m(x) m(x) • m(y)\ 
Vm(y) • m(x) m(y) • m(y)y 

( Gd(0) Gd(x - y)\ 

l^Gd(x-y) Gd(0) ) 


(40) 

(41) 


For d > 2 (the case d = 1 was studied in [5]), introducing 


In this case the correlation function is therefore 

C'r(x-y) =-^Gd(x-y). (35) 

Being the average optimal cost of the matching in 
the grid-Poisson case 


. Cd{G) 

detS(x,y)’ 

(42a) 

^ Cd(x-y) 

detS(x,y)’ 

(42b) 


_ ^(2) ^ (3g) 

we expect that in this case it will be asymptotically 
one half of the PP case. 

We will consider also the correlation function for the nor¬ 
malized transport field, i.e. the following quantity: 


observe that ^ 0+ for —>• oo, being NCd{x) finite 

for X 7 ^: 0 and NCd{0) ^ for d > 2, NCd{0) ~ In 

for d = 2. We have therefore that, in the notation above. 


det S = 


1 

A2 - B2 


(43) 


Cd(x-y) =cr(x) •(T(y), (37) and 


„ 1,, _ f A 2 'n -2 2ir 2 f j a „;„d—2 d D fO® j _ f°° .1 ™ _<i—l_d—1 4 cos 6> 

Cd(.x,yj — —J r(^-^ ) m d^sm ycos^Jp dmij^ dTO2TOi e ^ i 2/ 


B2r2 (^^) 

drHI) 




2 F 1 


52 


d+l d+1 



Cdi^^-y) 


(44) 


In the previous expression, we have introduced the hy¬ 


pergeometric function 


2 El 


a b 
c 


00 


-E 

n—O 


{a)n{b)n 
(c)„ n! ’ 


iS^')n 


r(Q-n) 

r (a — n -I- 1 )' 
(45) 
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Observe that we can reproduce exactly the same calcu¬ 
lation for the normalized field in the gP, obtaining 


r(i) 


E 


(2;gP) 

N 


id) 


Finally, for d > 2 we can compute also the so called wall- 
to-wall correlation function for the PP case: 


Wdir) := nf=2 Cd(r,X2, ■.■,Xd) 


jL I 

dN 


r(f) 


Gijr) 

Eiffd)' 


(47) 


Similarly, the computation for the gP case gives 
Wfir) := llt=2 cf(r,X2, ■■■,Xd) 


2 /r(^) 

dN 


r(i) 


Gi(r) 


(48) 


E 


(2;gP) 


(d) 


III. NUMERICAL RESULTS 

In the following we consider explicitly the cases d = 1, 
d = 2 and d = 3 and we numerically verify the results 
presented above. 


A. Case d — 1 


For d = 1 we have that 

I 'TTI 


Gi(r) = -E 


n^O 




1 |r| 

= -Y^ + y(l-H). (49) 


It follows from Eq. (32) that 


Ciix-y) = jf[^-\x- y\ (1 - |a: - y\)] ; (50) 

moreover, note that the average optimal cost is given by 

(51) 




In the gP case we obtain from Eq. (351 


Cfix-y) = 


(52) 


The average total cost of the optimal matching is given 
by 


E 


(2;gP) 

N 


(l) = Cf(0) = ^. 


(53) 


The previous results for the d = 1 case are known to 
the literature [DIHIIS], although the correlation function 
was derived using a different probabilistic approach. In 
[8] the correlation function for the normalized transport 
field is evaluated as 


Ciix) = Cj (x) = - arctan 


l-6a;(l-a;) 


12a:(l —a:)(l —3a:(l —a:)) 


( 54 ) 


B. Case d — 2 

For d = 2, denoting by x = (xi, X 2 ) and by y = (j/i, 2 / 2 ) 
two points on the unit flat torus T^, the Laplacian 
Green’s function can be written in terms of special func¬ 
tions as [H] 


G2(x-y) = -X; 


47^2 


1_^ p27rin (x-y) 


= 2F In 


o_3/4 i?i(^z|i) 

r(i) 




2=(xi-yi)+i{^2-y2) 


(55) 


where we introduced the first Jacobi theta function 

00 

di(z|T) := 2 e^ ^(-1)” e“™("+b gjn [(2n + l)z] 


n=0 


(56) 


From Eq. (32) we have simply 


1 

— In 

0 3/4^1 (ttiN) 

xl 

N 

27T 

r(f) 

2 

Z — Xi-\-iX2 

( 


In the gP case, we have as usual 


C 2 (x-y) = -C 2 (x-y). 


(58) 


Observe that the previous expressions contains no free 
parameters and therefore a direct comparison with nu¬ 
merical data is possible. We present our numerical re¬ 
sults both for the gP case and the PP case in fig. 

The average optimal cost for the PP Ebmp is given by 
e)^\2) = 02 ( 0 ): however, G 2 (x) is divergent for x = 0. 
Analysing the scaling of the total cost and performing a 
proper regularization of the previous quantity, Caracciolo 
et al. [ 1 ] obtained: 


G2(0) 


= 11 


In A 


N \ 2tt 


■dr 


+ o(^), /?PP = 0.1332(5). 


(59) 

A numerical fit of the optimal costs for d = 2 for the gP 
Ebmp gives 


cno) = 


1 /InA 


■ /3gp^ +0 ^ 


, /?<,P = 0.3758(5). 

(60) 

The correlation function (|37|) for the normalized 


matching field in the PP case has the expression ( |44[ ), 

7rG2(x-y) 

o4x-y) = ^ C 40 ) ■ 

Observe that the only free parameter in this quantity is 
^ 2 ( 0 ): inserting the value obtained by Caracciolo et al. 
[1], Eq. (591, we obtain the theoretical prediction in 
fig. |2b[ where we also present some numerical results for 
C 2 (x, y) that show the agreement with the theoretical 
curve. 
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a(p) P fixed 

p = 1 

0.860(5) 

to 

0.996(5) P = Ppp = 0.1332(5) 

p = 3 

0.96(1) 

p = 1 

0.88(1) 

gP p = 2 

1.02(3) P = Pov = 0.3758(5) 

p = 3 

0.98(2) 


Table I. Fit results for the wall-to-wall correlation function 
of the Euclidean matching problem for d = 2 and p = 1, 2,3 
with reference to the notation of the fitting curve, Eq. (64l. 
Observe that we expected q.{2) = 1. 


if we are dealing with a PP matching. We expected that 
a(2) = 1. We performed a numerical fit using the pre¬ 
vious expression also for p = 2, obtaining the results 
presented in Table Numerical data are in excellent 
agreement with the functional expression (64), suggest¬ 
ing therefore that the wall-to-wall correlation function in 
d = 2 has the same expression for all values of p up to 
a non universal multiplicative constant depending on the 
exponent that appears in the weight function. However, 
further investigations in this direction are needed to con¬ 
firm this results in a wider range of values of p. 


Using the value (60) in Eq. (46) for d = 2 we also 


obtained the theoretical curve for the grid-Poisson prob¬ 
lem depicted in fig. ( [2b| ), where, once again, an excellent 
agreement is found with numerical data. 

Finally, let us compute the wall-to-wall correlation 
function for the PP Euclidean matching problem on the 
flat torus. The theoretical prediction is given by Eq. (47), 


W2{r) = - 


2NC2{0) 


Gi(r). 


In the gP case, instead, we have 


(62) 


wr (r) = - 


4iVGfP(0) 


Gi(r). 


(63) 


Numerical results both for the PP case and for the gP 
case are presented in fig. and fig. |2d| Once again 
observe that the values of the average optimal cost in the 


corresponding cases, Eq. (^ and Eq. ( |59[ ), fix completely 
the expression of the wall 


-to-wall correlation function. 


The case p ^ 2 


Up to now, we analyzed the correlation functions for 
p = 2 and d = 2 obtained from the linearized equation 
(18). We present some numerical results for p = 1 and 


p = 3 and N = 3600 both in the gP case and in the PP 
in the two dimensional case. In particular, we analyzed 
the wall-to-wall correlation function for different values 
of p and we obtained numerical evidences of a functional 
form for it of the type 


W 2 {r\p,N) = w{p,N) [i - r(l - r)] 


(64) 


both in the PP and in the gP case. Inspired by the ob¬ 
tained expression for the p = 2 case, Eq. (62), we assumed 


for the global factor w{p) the following dependence on the 
size N of the considered system 


w{p,N) = 


»:(p)7 


2 In -|- 47r/3 


(65) 


C. Case d = 3 


The general expressions, Eq. (32) and Eq. (pfSl), for the 


correlation functions presented above can be applied in 
any dimensionality in the iV —>■ oo limit. We show here 
some numerical results for d = 3 for the grid-Poisson 
case, taking however into account that the precision of 
these data is lower, since the computation has complexity 
O(iV^) (i.e., the computational complexity of the Hun¬ 
garian algorithm in the Edmonds and Karp [4] version), 
where the number of points N has to scale exponentially 
with the system dimensionality in order to maintain the 
same accuracy. 


For d = 3 Eq. (22) and Eq. (32) give 


1 


1 


G3(x-y) = -^ y (66) 

' ' 2Tr^N ^ llnll^ 

neZ3\{0} 


Clearly the previous function can not be represented in 
a plot. From the correlation function G 3 (x), the wall to 
wall correlation function can be obtained as before in the 
form 


^3(0 = 


(67) 


As in the previous cases, G 3 ( 0 ) can be evaluated from the 
cost fit [I] and it is equal to G 3 ( 0 ) = 0.66251(2)A^“3 — 
(note that an exact formula for the coefficient 
of the ^ correction to G 3 ( 0 ) is provided in [T] in terms 
of an Epstein function). 

Following the same procedure of the PP case, we can 
compute the wall-to-wall correlation function on the unit 
hypercube in d = 3 for the gP matching problem. Re¬ 
producing the computations of the d = 2 case we have 

= ( 68 ) 

We evaluated G3^(0) from the cost scaling, obtaining 

GfP(O) = 0.4893(4)iV-i - (69) 


where a{p) depends only on the weight exponent p and The prediction obtained and the numerical data are pre- 

/? = /3gp if we are considering a gP matching, /? = /3pp sented in fig. 
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n 

(a) Section C2(ri,0) and C2'^(ri,0) of the correlation 
function both in the PP case for N = 10"* and in the gP 
case for N = 3600 and corresponding theoretical predictions. 



n 

(b) Section C2(ri,0) for N = 10^ and C2^(ri,0) for 
N = 3600 of the correlation function and theoretical 
predictions, Eq. \AA\ and Eq. ( |46[ |: note that the 
theoretical curves overlap. 



(c) Rescaled wall-to-wall correlation function in two 
dimensions for the PP matching problem with N = 3600 on 
the unit fiat torus. The continuous line corresponds to the 
analytical prediction. 



(d) Rescaled wall-to-wall correlation function in two 
dimensions for the gP matching problem with N = 3600 on 
the unit flat torus. The continuous line corresponds to the 
analytical prediction. 


Figure 2. Theoretical predictions for the correlation functions C 2 (x) and C 2 (x) for the Euclidean bipartite matching problem in 
two dimensions and numerical results. Numerical results for the wall-to-wall correlation functions and corresponding theoretical 
predictions are also presented. 


IV. CONCLUSIONS 

In the present work we adopted the scaling ansatz pro¬ 
posed by Caracciolo et al. [T] to compute the correlation 
function for the optimal matching ray and for the nor¬ 
malized optimal matching ray in the Euclidean bipartite 
matching problem on the d-dimensional flat hypertorus 
with quadratic cost. We showed also that the correlation 
function is strictly related to the Green’s function of the 
Laplacian operator on the flat hypertorus itself in the 
large N limit. Given the value the average optimal cost 
at fixed size N, the obtained expressions have no free pa¬ 
rameters and were directly compared with the results of 


numerical simulations, showing an excellent agreement. 
For d = 2 and d -|- 3 we computed also the wall-to-wall 
correlation function: for d = 2 in particular we give nu¬ 
merical evidences that, for p 2, the wall-to-wall corre¬ 
lation function has the same form obtained for p = 2 up 
to a global multiplicative constant. 

All previous results suggest that the Ebmp with 
quadratic cost, in the large N limit, appears as a Gaus¬ 
sian free theory on the d-dimensional flat hypertorus, in 
such a way that the correlation function of the matching 
ray is related directly to the free propagator of the theory 
itself. In subsequent publications we will investigate this 
crucial aspect of the problem and its implications on the 
















































Figure 3. Wall-to-wall correlation function in three dimen¬ 
sions for the gP matching problem with d = 3 and N — 9261 
on the unit flat hypertorus. 


universal behavior for different values of the exponent p 
in the cost functional. 
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